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The measurement of cosmic ray energy spectra, in particular for individual species, is an essential approach in finding 
their origin. Locating the “knees” of the spectra is an important part of the approach and has yet to be achieved. 
Here we report a measurement of the mixed Hydrogen and Helium spectrum using the combination of the ARGO-YBJ 
experiment and of a prototype Cherenkov telescope for the LHAASO experiment. A knee feature at 640 ± 87 TeV, 
with a clear steepening of the spectrum, is observed. This gives fundamental inputs to galactic cosmic ray acceleration 
models. 
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Introduction Galactic cosmic rays are believed to originate 
at astrophysical sources, such as supernova remnants. The 
mechanism for accelerating nuclei to energies from 10 14 eV 
to 10 2 ° eV remains unknown. A handful of significant struc¬ 
tures in the approximately power law spectrum occur over the 
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whole energy ranged. One of them is a significant down¬ 
ward bending of the spectrum around 3 x 10 15 eV, the so- 
called “knee”—. Many acceleration models have successfully 
explained the power-law characteristics of the spectrum, al¬ 
though no originating source has yet been experimentally ob¬ 
served for the high energy particles^. The knee of the spec¬ 
trum obviously plays a key role to test the proposed accel¬ 
eration and propagation models. One of the theories is that 
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the knee marks the highest energy that the galactic cosmic 
ray sources can reach-. The spectrum of all cosmic rays, 
however, does not appear to bend sharply, because different 
species may have different cut-off energies and extra-galactic 
cosmic rays may merge into the flux. These latter may dom¬ 
inate the flux at higher energies! Such a straightforward in¬ 
vestigation unfortunately has been very difficult in the past 
decades due to two experimental limitations. 1) the precision 
measurement of cosmic ray species and energies with space or 
balloon-borne calorimeters and charge-sensitive detectors has 
been constrained by their small exposure due to limited pay- 
load, so that statistically reliable measurements cannot effec¬ 
tively extend to energies higher than 10 14 eV—, far below the 
knee. 2) Ground-based experiments with extensive air shower 
(EAS) techniques are troubled by large uncertainties such as 
unknown energy scale and lack of effective tools to tag the na¬ 
ture of the primary inducing the observed showers, indepen¬ 
dently of the statistical accuracy of the measurement!!. As a 
consequence, different experiments find a different knee en¬ 
ergy as summarized in FIG. 1 of reference! mainly because of 
the unknown mixture. The uncertainty in the attempts of mea¬ 
suring the pure proton spectrum is still large, e.g. the knee is 
found a few hundreds of TeV in CAS A-MIA!! and a few PeV 
in KASCADE!. The lack of well-measured knee energies for 
individual species is prohibitive for developing a precise the¬ 
ory about the origin of cosmic rays. 

The situation is improved by the ARGO-YBJ experiment, 
at 4300 m above sea level in Tibet, which records nearly ev¬ 
ery secondary charged particle of showers incident upon its 
unique detector made of a continuous array of Resistive Plate 
Chambers (RPC)!!. Such a set-up brings the threshold of the 
shower measurement by Argo-YBJ down to the same energy 
range of CREAM!. This enables ARGO-YBJ establishing 
the energy scale by measuring the moon shadow!! and cross¬ 
checking with CREAM-M. This improvement is enhanced 
with the addition of data from a Cherenkov telescope- 14 imag¬ 
ing every shower in its field-of-view (FoV). The hybrid of 
the two techniques improves the resolution for shower energy 
measurements, and enhances the capability to discriminate 
showers induced by Hydrogen and Helium nuclei (//&//( ) 
from events initiated by heavier nuclei!!. Here, we report 
the measurement of the knee of the energy spectrum of the 
light component (HSzHe) below 1 PeV using the hybrid data 
from the ARGO-YBJ RPC array and the Cherenkov telescope, 
which is a prototype of one of the main instruments in the fu¬ 
ture LHAASO experiment^!!!. 

The Hybrid Experiment The hybrid experimental data set 
includes air showers whose cores are fallen inside an area of 
76 to x 72 m fully covered by the ARGO-YBJ RPC array, i.e., 
1 m from the edges of the array, and whose arrival directions 
are in the effective FoV of the telescope, i.e., a cone of 6° 
with respect to the main axis of the telescope, which has a full 
FoV of 14° x 16° pointing to 30° from the zenith. The tele¬ 
scope is about 79 to off the center of the array in the south-east 
direction. This defines a geometrical aperture of 163 m 2 sr. 
According to the simulation of the hybrid experiment, high 
energy (> 100 TeV) showers are detected with almost full ef¬ 
ficiency, particularly the HSzHe events. This minimizes the 


uncertainty of the cosmic ray flux measurement. 

In its FoV, the telescope- has an array of 256 pixels with 
a size of approximately one square degree each. The shower 
image records the accumulated Cherenkov photons produced 
in the entire shower development. As described below, the 
total number of photons in the image can be used to recon¬ 
struct the shower energy. The image shape indicates the depth 
of the shower development after reaching its maximum, giv¬ 
ing useful information to select proton or Helium showers. 
The ARGO-YBJ array consists of 1836 RPCs, each equipped 
with two analog readout “Big Pads” (140cm x 123cm) to col¬ 
lect the total charge induced by particles passing through the 
chambe r 1 *- 1 *' 1 . The collected charge is calibrated to be propor¬ 
tional to the number of charged particles^—. The most hit 
RPC, together with the surrounding RPCs, measures the lat¬ 
eral distribution of secondary particles within 5 m from the 
core. Such a unique measurement is very useful not only for 
a precise reconstruction of the shower geometry, but also for 
the selection of proton and Helium showers. 

The coincident cosmic ray data, collected in the hybrid ex¬ 
periment from December 2010 to February 2012, are used for 
the analysis presented in this paper. The main constraint on 
the exposure of the hybrid experiment is the weather condi¬ 
tion in the moon-less nights. The weather is monitored by us¬ 
ing the bright stars in the FoV of the telescope and an infrared 
camera covering the whole sky. More details about the criteria 
for a good weather can be find elsewhere!!!!. Combining the 
good weather conditions and the live time of the RPC array, 
the total exposure time is 7.28 xlO 5 seconds for the hybrid 
measurement. Further criteria (quality cuts) for well recon¬ 
structed showers in the aperture of the hybrid experiment are 
1) at least 1000 particles recorded by the ARGO-YBJ digital 
readout)! to guarantee high quality reconstruction for shower 
geometry*!; 2) at least 6 pixels triggered in each shower im¬ 
age. About 32,700 events survived these cuts. Among them, 
8218 high energy events approximately above 100 TeV are 
detected. The core and angular resolutions are better than 1.2 
to and 0.3°, respectively. 

A great deal of extensive air showers, including their 
Cherenkov photons, are simulated by using the CORSIKA 
codec! with the high energy hadronic interaction model 
QGSJETII-03*! and with the low energy model GHEISHA*!. 
The G4argo* 7 package and a ray-tracing procedure on the 
Cherenkov photons*! are applied for further simulation of the 
detector responses. All five mass groups, i.e. proton, Helium, 
C-N-O group, Mg-Al-Si group and Iron are generated in the 
simulation. A detailed comparison between the data and the 
simulation can be found elsewhere!!. 

Shower Energy Reconstruction and All-Particle Distri¬ 
bution The shower energy, E, is reconstructed using the total 
number of Cherenkov photo-electrons, N pe , collected by the 
telescope which observes the shower with a certain impact pa¬ 
rameter R p , given by the shower geometry. Using a very large 
sample generated by the simulation described above, a look¬ 
up table for the shower energy with two entries, i.e., N pe and 
R p , is determined for each mass group. For a shower with 
N pe measured by the telescope and R p measured by the RPC 
array, the shower energy can be read out from this table. The 
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FIG. 1. Distribution of the number of Cherenkov photo-electrons 
measured by the telescope (filled circles) with a bin width of 0.22 
in logioN pe . The histograms represent the distributions generated 
according to the flux models- 30,31 and to the all-particle spectra and 
corresponding composition models obtained from Tibet AS-y— and 
KASCADB^ experiments. 


energy resolution is found symmetric and fits well a Gaus¬ 
sian function with a between 23% and 27% for different mass 
group. However, a clear feature of the energy reconstruction 
is a systematic shift which depends on the nature of the pri¬ 
mary. Around 1 PeV, the difference between proton and Iron 
showers is approximately 37%, significantly greater than the 
resolution. For a mixed sample with unknown composition, 
this feature will distort the all-particle energy spectrum even 
if the measurement is fully efficient, namely with a constant 
geometrical aperture independent of the energy. 

In order to compare with other experiments or existing cos¬ 
mic ray flux models without assuming any specific mixture 
of the species for the measured showers, in FIG. Q] we plot 
the event distribution as a function of the measured number 
of photo-electrons N pe . Also plotted in the same figure are 
the distributions generated according to the all-particle spec¬ 
tra measured by the experiments and the corresponding as¬ 
sumptions on the mixture of different species. Here we show 
the results from Tibet AS 7 with two different composition 
models^, from KASCADE with its composition models ob¬ 
tained from the unfolding procedures^, and two widely quoted 
composition models, i.e., Hdrandel—! and H4a— The corre¬ 
sponding energy range is from 126 TeV to 15.8 PeV assum¬ 
ing 1:1 mixture of proton and Helium showers. The compar¬ 
ison shows that the existing all-particle spectra and their cor¬ 
responding composition models are in a general agreement at 
a level of 30%. The data used in this work also maintains a 
general agreement with others at a similar level. 

Hydrogen and Helium Event Selection The secondary 
particles in showers induced by heavy nuclei are spread fur¬ 
ther away from the core region. Therefore significant differ¬ 
ences of the lateral distributions exist in the vicinity of the 
cores between showers induced by light or heavy nuclei? 1 . 
Beyond a certain distance, e.g., 20 to from the core, the 
lateral distributions become similar because they are mainly 
due to multiple Coulomb scattering of the secondary parti¬ 
cles and are well described by the Nishimura-Kamata-Greisen 
(NKG) function. With its full coverage, the ARGO-YBJ ar¬ 


ray uniquely measures the lateral distribution of the secondary 
particle density at the shower core. The number of particles 
recorded by the most hit RPC in an event, denoted as N rnax , 
is a good parameter to discriminate between showers with dif¬ 
ferent lateral distribution within 3 m from the cores. In a 
shower induced by a heavy nucleus, N max is expected to be 
smaller than that in a shower induced by a light nucleus with 
the same energy. Obviously, N max depends on the shower 
size which can be indicated by N pe at the distance R p from 
the shower axis. We define a reduced dimensionless variable 
Pl = log w N max - lA4log w N pe - R p /81.3m + 3.3 empiri¬ 
cally obtained from the MC simulation^ to absorb the shower 
size effect. 

The shape of the shower image recorded by the Cherenkov 
telescope is also a mass-sensitive parameter. The elliptical 
image is described by the Hillas parameters^, such as width 
and length. The images are more stretched, i.e., narrower 
and longer, for showers that are more deeply developed in 
the atmosphere. The length to width ratio ( L/W ) is there¬ 
fore a parameter sensitive to the depth of the shower maxi¬ 
mum which depends on the nature of the primary. It is also 
known that the images are more elongated for showers farther 
away from the telescope, because of purely geometric rea¬ 
sons. The ratio L/W is nearly proportional to the shower 
impact parameter R p , but depends very moderately on the 
shower size. Taking into account the dependence on mea¬ 
sured number of Cherenkov photons in a shower and on the 
impact parameter, we define a reduced dimensionless variable 
pc = L/W — R p /97.2m — 0A4log\aN pe + 0.32, obtained 
again from the MC simulation-^, to absorb both the R p and 
shower size effects. 

The selection of the II&IIe sample is carried out by com¬ 
bining the two composition-sensitive parameters. A correla¬ 
tion analysis shows that the two variables are quite indepen¬ 
dent. Composition groups are significantly separated in the 
Pl-Pc rnapii. This map can be plotted with probability con¬ 
tours for the two mass groups, namely H&He and the other 
nuclei (FIG. 0. The cuts p l > —1.23 or pc > 1.1 result 
in a selected sample of //&//e showers with a purity of 93% 
below 700 TeV and an efficiency of 72% assuming the compo¬ 
sition models given i n 30 ' 31 . The aperture gradually increases 
to 120 m 2 sr at 300 TeV and keeps as a constant at higher en¬ 
ergies. The contamination of heavier nuclei increases with en¬ 
ergy, becoming about 13% around 1 PeV and gradually rising 
to 45% around 6.5 PeV. The contamination obviously depends 
on the composition assumption^. The associated uncertainty 
is discussed below. 

Energy Spectrum of Proton and Helium The energy 
spectrum of the selected sample of H&cHe showers is plotted 
in FIG. [3] The shower energy is now better defined because 
the intrinsic scale difference between H and He showers is 
smaller than 10%, significantly lower than the energy resolu¬ 
tion. Using about 40,000 simulated events that survived all 
the reconstruction quality cuts and TT&fJe selection, the en¬ 
ergy resolution function is found to be Gaussian with a con¬ 
stant standard deviation of 25% and overall systematic shift 
less than 2% at energies above 100 TeV—. To take into ac¬ 
count the energy resolution and any kind of smearing like bin- 
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FIG. 2. Composition-sensitive parameters pl and pc for the two 
mass groups, H&zHe (solid contours) and heavier components 
(dashed contours). The numbers on the contour isolines indicate the 
percentage of contained events. 


to-bin migration from the true to the reconstructed primary 
energy, a Bayesian algorithm^ is applied to unfold the ob¬ 
servational data. The selection efficiency for He showers is 
approximately 80% that for H showers. The contamination 
due to heavy nuclei is subtracted in each bin considering the 
composition model in Ref. 30 . 

Systematic uncertainties mainly arise from the following 
causes: 1) Assumptions on the flux of heavy species. Below 
800 TeV, a flux uncertainty of 1.9% is estimated by consid¬ 
ering the models in ref. 30 , ref^i and the extrapolation from 
the CREAM data^. Above 800 TeV, this uncertainty increases 
with energy up to 23.3% at 2.5 PeV. 2) Because of slightly dif¬ 
ferent detection efficiencies for H and He showers, the frac¬ 
tion of Helium in the selected samples depends on the com¬ 
position assumption. This also results in an uncertainty of 3% 
in the overall flux. 3) Choice of the interaction models. The 
overall flux uncertainty is about 4.2% by considering the high 
energy interaction models SIBYLLA and QGSJET, and the 
low energy interaction models GHEISHA and FLUKA—f 4) 
Boundary effect of the aperture. The boundaries were slightly 
varied to smaller RPC array and smaller FoV of the telescope. 
The corresponding flux uncertainty is about 3% indicated by 
the deviation from a linear response to the variation. 5) Cal¬ 
ibration for the number of particles measured by the RPCs. 
Depending on the calibration methods, we find an uncertainty 
of 7% in the number of events surviving all the criteria which 
involve the RPC response. The overall systematic uncertainty 
on the flux is plotted as the shaded area in FIG. [3] 

A more strict cut for higher purity (97%) HSzHe sample 
has been applied below 700 TeV where the spectrum fits well 
with a single-index power law, according to CREAM- and 
ARGO-YB.U-. This yields a much smaller but constant aper¬ 
ture of ~ 50m 2 sr above 250 TeV, a negligible contamination 
from heavy nuclei and a corresponding precise measurement 
of the spectrum^ shown by the filled squares in FIG. [3] This 
serves as a verification for both energy scale and absolute flux 
once compared with the previous measurements. The differ¬ 
ences between the fluxes measured by the above experiments 


is found less than 9%H. 

Discussion and conclusions An evident bending struc¬ 
ture is observed in the cosmic // & H e spectrum by the hy¬ 
brid experiment using the ARGO-YBJ RPC array and the 
Cherenkov telescope at 4300 m above sea level. The previ¬ 
ous measurement s 5 ' 13 ’ 15 below 700 TeV, as mentioned above, 
indicate that the spectrum is a single power law with index 
-2.62. Beyond 3 PeV, however, many experiments reported 
an evident bending in the all-particle spectrum, although with 
a spread in the bending energy. With the normalization fac¬ 
tor Jo = (1.82 ± 0.16) x 10 -11 GeV~ x m~ 2 s _1 sr _1 at 
400 TeV the four flux values measured below 700 TeV are 
fitted well by the single-index power law. Together with its 
extrapolation up to 3160 TeV, the upper boundary of the last 
bin in our analysis, the power law spectrum has been used as 
the a priori expectation. The significance of the deviation of 
the bent spectrum from the single-index power law is mea¬ 
sured by calculating the chance probability as follows. 85 
events in the range from 800 TeV to 2000 TeV and 9 events 
from 2000 TeV to 3160 TeV are observed, to be compared 
with an expectation of 118 and 22 events, respectively, de¬ 
rived from the hypothetical spectrum plus the contaminating 
heavy species, corresponding to a deviation of 4.4 a. A bro¬ 
ken power law fits well the measured spectrum. Below a break 
energy Ek, the assumed spectrum describes the data very well 
with x 2 /dof = 0.7 for the first four points. Above the break 
the data can be fitted by Jq ■ (f?fc/400TeV) _2 62 • ( E/Ek ) /3 , 
with E k = 640±87 TeV and /3 = -3.34±0.28. The relatively 
large error on E k is due to the limited statistics and the finite 
energy resolution. In addition, the systematic error in the en¬ 
ergy scale is 9.7%ii, which corresponds to ~ 62 TeV at /%— . 

In summary,the joint operation of the ARGO-YBJ detector 
with a wide field-of-view and imaging Cherenkov telescope 
allowed a detailed investigation of the energy range bridg¬ 
ing the gap between the direct observations of CREAM and 
the ground-based KASCADE experiment. This hybrid exper¬ 
iment yields a clear evidence for a knee-like structure in the 
spectrum of light primaries (protons and Helium nuclei) be¬ 
low 1 PeV.The observation of the knee of the primary light 
component at such a low energy gives fundamental inputs to 
galactic cosmic ray acceleration models. 
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